Shot Noise of a Quantum Shuttle 
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We formulate a theory for shot noise in quantum nanoelectromechanical systems. As a specific 
example, the theory is applied to a quantum shuttle, and the zero-frequency noise, measured by 
the Fano factor F, is computed. F reaches very low values {F ~ 10~^) in the shuttling regime 
even in the quantum limit, confirming that shuttling is universally a low noise phenomenon. In 
approaching the semiclassical limit, the Fano factor shows a giant enhancement [F ~ 10^) at the 
shuttling threshold, consistent with predictions based on phase-space representations of the density 
matrix. 
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Nanoelectromecanical systems (NEMS) are presently a 
topic of intense research activity [l| . These devices com- 
bine electronic and mechanical degrees of freedom to dis- 
play new physical phenomena, and potentially may lead 
to new functionalities. An archetypal example of such a 
new phenomenon is the charge shuttling transition exhib- 
ited by the device proposed by Gorelik et al. 0: here a 
movable nanoscopic object begins to transport electrons 
one-by-one beyond a certain threshold bias. Recent work 
has extended the original ideas to the quantum regime 
(the motion of the movable part is also quantized), and 
has shown that the shuttling transition occurs even in 
this limit, albeit in a smeared-out form 

An unequivocal experimental observation of the shut- 
tling transition has not yet been achieved. The IV-curve 
measured in the recent experiments on a Cgo single- 
electron transistor can be interpreted in terms of shut- 
tling [3, but also alternative explanations have been pro- 
motedj^, 01 ^3 ■ It is therefore natural to look for more 
refined experimental tools than just the average current 
through the device. An obvious candidate is the cur- 
rent noise spectrum [ill fl^ . The measurement of the 
noise spectrum or even higher moments (full counting 
statistics) reveals more information about the transport 
through the device than just the mean current. The the- 
oretical studies of the noise have attracted much atten- 
tion recently in NEMS in gene ral [isl lisj as well as 
for the shuttling setup [la IITI in the (semi)classical 
limit. More specifically, Pistolesi 16] reports a vanishing 
Fano factor in the large amplitude limit of a driven shut- 
tle. On the other hand, a study by Isacsson and Nord 
[itI i of a classical system exhibiting shuttling instability 
found, somewhat surprisingly, a higher Fano factor in the 
shuttling regime than in the tunneling regime. This re- 
sult is attributed to a different confining potential used 
in ^3 compared to the other studies in this area. 

The aim of this Letter is to develop a quantum me- 
chanical theory for the shot noise spectrum for quantum 



NEMS, and apply it to the model introduced by Gorelik 
et al. [21 • Our method combines the classical nature of the 
charge transfer processes in the high bias limit with 
an operator version of a generating function technique 
[2(tI |. While the present work considers only Markovian 
master equations, we believe that the method can be gen- 
eralized to the case where the dynamics of the mechanical 
degrees of freedom is non-Markovian. For systems where 
the current noise can be expressed in terms of system 
operators (using the quantum optics language), such as 
the quantum dot array of Ref. an alternative eval- 
uation of noise, based on the quantum regression theo- 
rem (QRT) is possible, and we have verified that the two 
methods give identical results in this case . We stress 
that the converse is not true: QRT is not applicable to 
the single-dot case. 

Our previous quantum calculation Q of the mean cur- 
rent relied on a generalized master equation (GME) for 
the system density matrices pii{t) (pn and poo describe 
the occupied and empty dot, respectively, and the off- 
diagonal components decouple from their dynamics, and 
can be neglected). In order to compute the noise spec- 
trum, we follow the ideas of Gurvitz and Prager and 
introduce number-resolved density-matrices pf^\ where 
n = 0, 1, . . . is the number of electrons tunneled into the 
right lead by time t. Obviously, pii{t) = Pu \'^)- "^^^ 
Pu^ obey 

{t) = ^ [^osc, p[,o^ (t)] + /:da,„p it) 

.E±{e-^^/\p(^)(^t)} + Tne^/^p^ri-'\t)e^^\ 

p[f{t) = l[i7o3c - eEx, p[l\t)] + /Idamp P^n\t) 

(1) 

with p[i^\t) = 0. In the commutators describe co- 
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herent evolution of discharged or charged harmonic os- 
cillator of mass m and frequency uj in electric field re- 
spectively. The terms involving F^ jj describe the charge 
transfer processes from/to leads while the mechanical 
damping with the damping coefficient 7 is determined 
by the kernel (at T = 0) H 
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[x,[x,p]]. (2) 



The mean current and the zero-frequency shot noise 
are given by 



/ = e-VnF„(t) =eVnP„W 
at ^ — ' ^ — ' 



5(0) = 2e^- 



(3) 
(4) 



where P„(t) = Trosc[Poo^(*) + are the prob- 

abilities of finding n electrons in the right lead by 
time t. Using Eq. we find / = J2n'"-^nit) = 
TiiTTosc{e^^^^piiit)), i.e. one recovers the station- 
ary current used previously 0. In a similar fashion, 

En^'Pnit) = r^Trosc{e2^/^[2E„np(';^(i)+Pii(0]}, 
whose large-time asymptotics determines the shot noise 



according to This can be computed using an 

operator-valued generalization of the generating func- 
tion concept. We introduce the generating functions 
Fii{t;z) = Y.nPu\^)^'^ "^ith the properties Fii{t]l) = 
Pii{t), ■§^_Fii{t] z)\:,=i = Y.n'^Pvi\^)- The equations of 
motion for Fii{t; z) are 



^Foo{t;z) = ^[H,,c,Fooit;z)] - lL{e-^-/\ Foo{t; z)} 
at in 2 

+ /:damp -Foo(i; z) + zTne^^^Fn {t; z)e"/^ 

= CooFoo{t; z) + zCQiFii{t;z), 

d 1 

— Fii{t;z) = — [Hose - eEx,Fii{t;z)] + Cda.mp Fii{t; z) 
at in 

-^{^^/\F,,{t-z)} 

+ TLe--^^Foo{t;z)e~^^^ 
= CioFQo{t; z) -f CiiFii{t] z), 

(5) 



where we have introduced the block structure of the Li- 
ouvillean (super)operator C = ( ^™ ) . Using the F's 
the shot noise formula can be rewritten as 



J 



= {Tro,e[e2^/'( 



e^-/^ 2— Fn(i;z) +Fn(i;l) 
* oz z=l 



— 2Tr 
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-/^i^n(i;l)lTroscf|- Vi^.^(i;z) 1} . (6) 
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A Laplace transform of jSJ yields 



Fm{s\ z) 
Fnis;z) 



s - Cqo -zCqi 
—CiQ s — L\\ 



where /j'f *(-z) = /^if' (0)2^" depends on the initial con- 
ditions. Defining the resolvent Qi^s) — ~ C)~^ of the 
full Liouvillean we arrive at 



(8) 



is singular at s = 0. To extract the singular behavior we 
introduce the projector V on the null space of the Liou- 



villean: V* 
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Trsys(»). We also need the com- 



plement Q = \ —V. Using the relations VL = CP ~ 
and L = Q£Q, the resolvent can be expressed as Q{s) = 
{sP + sQ~QCQ)-^ = IV+Qj^Q ^ \P-QC-^Q, in 
leading order for small s. The object QC~^Q (the pseu- 
doinverse of C) is regular as the "inverse" is performed on 
the Liouville subspace spanned by Q where C is regular 
(no null vectors). 
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In order to extract the large-i behavior wc study the 
asymptotics of the above expressions as s ^ 0+. This is 
entirely determined by the resolvent Q{s) in the small-s 
limit. Since C is singular (recall Cp^*"^^ = 0) the resolvent 



Substituting the asymptotic behavior of the resolvent 
into Eqs. © and (|5J|, keeping only the terms divergent 
at s = in both equations, and performing the in- 
verse Laplace transform |25| we find the following large-t 



asymptotics 
d_ fFooit;z) 
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where we have defined an auxihary quantity S — 
Q/:"^Q(r«'=°'^^Pi'l°''="'^) and C™* depends on initial 
conditions. Using these in ((HJ we arrive at the final ex- 
pression for the Fano factor F = S{0)/2eI 

^^ = l_MiiTro,e(e2^/'Sii). (11) 

It is of crucial importance that this expression is inde- 
pendent of the initial conditions [in the algebra leading 
to Hll(l the linearly divergent terms in t and the initial 
condition terms cancel identically] . S satisfies 



"P ^x/ A „stat 'X /A / „stat 



/ „stat 

:Pii 



with TrsysS = 0. 



(12) 

Equations l|ll|) , p2|l together with the stationary version 
of the GME 



^(^|It ) =0, with Tr,y,p^*^* = 1 



(13) 



form the main theoretical result of this Letter and are the 
starting point for the calculation of the noise properties 
of the quantum shuttle. 

In general, these equations have to be solved numer- 
ically. However, there is an analytic solution to them 
in the limit of small bare injection rates compared to 
damping, i.e. ^ ^ 7 <C w. In this limit the os- 
cillator gets equilibrated between rare tunneling events 
and, consequently, the state of the oscillator in a given 
charge state is close to its corresponding canonical state, 
i.e. pfo^'^ = PooPosc(O), pff* = piiPosc(e-E), posc(0 = 
e-'^(^°-^'^VTrosce-^(^°--''^) where only the probabil- 
ities Poo, Pi 1 of the respective occupations are to be de- 
termined from (|13|l . By tracing Eq. H13|) with respect 
to the oscillator we find — F^poo + ^RPii = (with 
Poo+Pii = 1) where = riTrosc(e-2-/Vo«c(0)), = 
rijTrosc(e^^/^Posc(ei?)) are the renormalized tunneling 
rates. Proceeding similarly in the evaluation of E (this 
intuitive approach can be easily justified with singular 
perturbation theory, see e.g. Ref. |13l)) we set So*J** — 
sooPosc(O), Sff* = siiPosc{eE) and tracing Eq. ^ with 
respect to the oscillator we arrive at — F^soo + Fijsn = 
-Pii (with sqo + sii — 0). Solving the equations for 



Pa, Sii we find I = e 



F = 



which is 



the standard result for the two-state sequential tunnel- 
ing process . The Fano factor only depends on the 
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FIG. 1: Current / (upper panel) and Fano factor F (lower 
panel; log scale) vs. damping 7 for different transfer rates F 
and tunneling lengths A. The parameters are A = 2:0, F — 
O.OltJ (dot-dash-dashed line); A = a:o,r — 0.05tJ (dot- 
dashed line); A = 2xo,r = O.Oluj (full line); A = 2xo,V = 
0.05tj (dashes) with xo = ^yh/muj. Other parameters are 
eE/muj^ = 0.5a;o, T — 0. The asterisk defines the parameters 
of Wigner distributions in Fig. |21 



ratio between the rates: ^ = ^-'^P \m^'^ ■ This result 
provides us with the analytic asymptotical expressions for 
the current and the Fano factor for small hopping rates 
^L.R which we used to check our numerical routine. 

In the general case Eqs. (I13f) and (|12|l must be solved 
by truncation of the oscillator Hilbert space by retain- 
ing the N lowest states of the oscillator, and solving nu- 
merically the resulting 2N'^ x 2iV^ linear systems. Since 
the required N for a satisfactory convergence could reach 
N = 100 resulting in big (non-sparse) linear systems we 
employed the Arnoldi iteration and generalized minimum 
residual method (GMRes) for the solution of Eq. H13|) 
and (|12|l . respectively [i^. When properly implemented 
these methods provided a fast solution with a modest 
memory requirement j2ll |. 

In Fig. n we present the plots of the mean current and 
the Fano factor as functions of the damping coefficient 
for different parameters A, and F = F/, = F/{. The / 
vs. 7 curve shows the tunneling to shuttling crossover as 
damping is decreased, as explained in our previous work 
4]. The crossover spans a narrower range of 7's in case 
of A = 2xo {xq ~ ^hjraLo) compared to the A = a;o 
case. Thus, already for A = 2x0 the shuttle behaves al- 
most semiclassically, where a relatively sharp transition 
between the two regimes is expected. Around the transi- 
tion the tunneling and shuttling regimes may coexist, as 
shown analytically in We see this phenomenon ex- 



4 




FIG. 2: Phase space picture of the shuttle around the tran- 
sition where the shuttling and tunneling regimes coexist. 
From left to right we show the Wigner distribution func- 
tions for the discharged (Woo), charged (Wn), and both 
(Wtot = Wdo + Wii) states of the oscillator in the phase space 
(horizontal axis - coordinate in units of xo — ^/h/mu, vertical 
axis - momentum in h/xo). The values of the parameters are 
(corresponding to the asterisk in Fig.0: A — 2x0, eE /mu)'^ = 
0.5x0,7 = 0.029a;, r = Q.QIlo,T = 0. 

plicitly in Fig. [21 where we plot the Wigner distribution 
functions defined by 

(14) 

for a specific set of parameters corresponding to the 
"most classical" curve around the transition in Fig. 2] 
denoted by the asterisk. The Wigner plots show the 
quasiprobability distributions in the phase space of the 
oscillator resolved with respect to its charge state and 
prove the coexistence of the tunneling regime (character- 
ized by the spots around the phase-space origin) and the 
shuttling regime with the half-moon or ring-like shapes 
in Woo, Wii or Wtot, respectively fl- This semiclassical 
transition is accompanied by the nearly singular behavior 
of the Fano factor reaching the value « 600 at the peak. 
This is in agreement with the recent classical study 0| 
where the singularity of the Fano factor at the transition 
was also predicted. 

More important, however, is the behavior in the shut- 
tling regime. We can see in Fig. that the Fano factor is 
very small in the shuttling regime. This is true even in 
the strongly quantum case A = where the transition 
peak characteristic of the classical case is almost totally 
missing. As found previously j3| the classical transition 
is strongly smeared by the quantum noise into a broad 
crossover which is reflected by the absence of the peak in 
the Fano factor. Nevertheless, the shuttling regime still 
persists and is again characterized by a low noise. 

To conclude, we have presented a generic method of 
the calculation of the shot noise for quantum nanoelec- 
tromechanical systems and applied it to a quantum shut- 
tle system. We show that even in the quantum case 
the shuttling regime is characterized by a highly ordered 
charge transfer mechanism accompanied by the low cur- 
rent noise compared to the tunneling regime. When ap- 
proaching the semiclassical limit, the Fano factor shows 



a giant enhancement at the shuttling threshold, consis- 
tent with other classical studies and with the phase space 
analysis of the stationary density matrix. 
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